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Observation of the decay B° — > A+pri 
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In a sample of 467 million BB pairs collected with the BABAR detector at the PEP-II collider 
at SLAC we have observed the decay B — > A^pn and measured the branching fraction to be 
(1.94 ± 0.17 ± 0.14 ± 0.50) x 10~ 4 , where the uncertainties are statistical, systematic, and the 
uncertainty on the A£ — > pK~iv + branching fraction, respectively. We determine an upper limit of 
1.5 x 10~ 6 at 90% C.L. for the product branching fraction B(B° -> T+(2455)p) x B(A+ ->■ pK~ir + ). 
Furthermore, we observe an enhancement at the threshold of the invariant mass of the baryon- 
antibaryon pair. 

PACS numbers: 13.25.Hw, 13.60.Rj 



Although approximately 7% of B-meson decays have 
baryons in the final state, presently the sum of all mea- 
sured branching fractions of exclusive baryonic B decays 
is only about 1% [1]. B mesons decay dominantly via 
b — > c transitions, hence decays to baryons should be 
dominated by charm baryon production or a charmed me- 
son accompanied by non-charmed baryons. Both types 
of decays have been observed [2, 3], and are found to have 
comparable branching fractions for decays to final states 
with the same multiplicity. 

In baryonic B decays and in baryon production in 
general, enhancements at the threshold for the baryon- 
antibaryon invariant mass have been observed [3, 4] . This 
may indicate resonances near threshold or another mech- 
anism for enhanced production of baryon-antibaryon 
pairs. This threshold enhancement may also explain the 
increase in branching fraction with final state multiplic- 
ity and the apparent suppression of two-body decays to 
baryons [1, 5]. 

The mechanisms of baryon production in heavy meson 
decays are poorly understood, and studies of exclusive 
decays may provide insight into different decay mecha- 
nisms. As will be discussed below, isospin relations will 
also help distinguish different primary processes. 

In this paper, we present a study of the decay B° —$■ 
A+pir [14] and measure its branching fraction. The 
CLEO collaboration previously set an upper limit of 
B(B° — ► A+pn ) < 5.9 x 10~ 4 based on an integrated 
luminosity of 2.39 fb _1 [6]. For the isospin-related de- 
cay, B~ — > A+pTT~ , several measurements of the branch- 
ing fraction have been performed [7, 8]. The recent 
BABAR measurement gives (3.38±0.12±0.12±0.88)xI(T 4 
[9], a value that is significantly higher than earlier mea- 
surements (4.3cr deviation). The last and dominant error 
is due to the uncertainty in the A^ — > pK~ir + branching 
fraction, common to all measurements. 

While the B~ — > A+p-K~ final state can only have an 
isospin / of 3/2, B° ->■ A+pn° can also have / = 1/2. If 
both decays proceed via the same weak decay mechanism, 
/ = 3/2, the ratio of the partial decay widths of B° to 
B~ should be 2/3. However, it is also possible that the 
decay mechanisms are different. Thus a deviation of the 
ratio of partial decay widths from 2/3 would suggest a 
contribution from the / = 1/2 final state to the B° — > 
A~^pir° decay or a contribution from the decay process 
where the 7r~ is coming from the W in the B~ — > A^"pTT~ 
decay. 

This analysis is based on a dataset of about 426 fb~ 
corresponding to 467 million BB pairs. These data 
were collected with the BABAR detector at the PEP-II 
asymmetric-energy e + e~ collider with a center-of-mass 
energy, y/s, at the T(4S*) resonance mass. An additional 
sample of 44.5 fb _1 , collected 40MeV below the mass of 
the T(4S) resonance, are used to study the continuum 
background e + e~ — > qq, where q = u, d, s or c. 

The signal efficiency is determined using a detailed 



Geant4 [10] Monte Carlo (MC) simulation of the BABAR 
detector that generates MC events uniformly in the 
A+pir° phase space. MC events are also used to study 
the background contributions. 

The BABAR detector is described in detail else- 
where [11]. Charged particles are distinguished and their 
momenta measured in the tracking system consisting of 
a five-layer double-sided silicon vertex tracker (SVT) and 
a 40-laycr drift chamber (DCH). An internally reflecting 
ring imaging Cherenkov detector (DIRC) is also used to 
distinguish charged particles and a CsI(Tl) electromag- 
netic calorimeter (EMC) is used to detect photons. 

Likelihood ratios based on information from SVT, 
DCH and DIRC are used to identify protons and kaons. 
The efficiency for the kaon selection is around 90% while 
the rate for misidentifying pions and protons as kaons 
varies between 5% and 10%, depending on track momen- 
tum. The identification efficiency for the proton selec- 
tion is greater than 90% while the misidentification rate 
of identifying kaons and pions as protons varies between 
3% and 15%, depending on track momentum. 

Two photons are selected as electromagnetic showers 
in the EMC with the expected shape and are combined 
to form a 7r° candidate where the photon with the lower 
energy must have an energy greater than 60 MeV, while 
the second photon must have an energy greater than 
100 MeV. The invariant mass of the 77 combination 
is required to be between 120MeV/c 2 and 145 MeV/ c 2 . 
The yl+ candidates are reconstructed in the decay mode 
A£ — > pK~ir + , and a fit with geometric constraint ap- 
plied to the common vertex must have a x 2 probability 
greater than 0.1%. The invariant pK~n + mass must be 
within 2.5cr of the fitted peak of the mass distribution, 
2.276 < m(pK-n+) < 2.296 GeV/c 2 . The A+ and tt 
candidates are then combined with a p candidate in a fit 
using kinematic constraints to form a B° candidate. In 
the fit the mass of the pK~ir + candidate is constrained 
to the mass of the A+ and the mass of the 77 combina- 
tion to the mass of the ir° [1] . The x 2 probability of this 
fit must be greater than 0.1%. 

The analysis makes use of two almost independent 
kinematic variables, AE and Toes, where AE = E* B — 
y/s/2 is the difference of the reconstructed energy Eg and 
half of y/s in the e + e~ center of mass frame (CMS). The 
other variable is toes = \/(s/2 + p • Pb) 2 /-^ 2 — Pb 
where (Eg, po) is the four momentum of the e + e~ system 
and pb is the B candidate momentum, both measured 
in the laboratory frame. The toes distribution for sig- 
nal events peaks at the B mass and the distribution of 
AE for signal events is centered around zero. Candidates 
arising from other B decays, with more final-state par- 
ticles, such as B° — > A+pir + ir~ , are shifted to negative 
values of AE. Conversely, candidates arising from B de- 
cays with fewer final-state particles, such as B° — *■ A+p, 
are shifted to positive values. To suppress these decays 
only candidates with -50 MeV < AE < 40 MeV are se- 




FIG. 1: AE distribution for data signal events after all se- 
lection cuts (data points) and signal MC events (histogranr) 
normalized to the number of data signal events; signal events 
are obtained from binwise mES fits; dashed lines show the 
range used for tties distributions. 



lcctcd. 

A considerable background comes from B~ — > A~^pir~ 
decays, and in particular from the B~ — } Z , °(2455)/j, 
i7°(2455) — > A+ir~ decays, in which the A+p pair from 
B~ decay is combined with a tt° from the decay of the 
B + meson. To suppress this background, we reconstruct 
B~ -» A+pn" and reject the event if \AE\ < 50MeV 
and 771es > 5.27GeV/c 2 for such a B~ candidate or if 
the condition 2400MeV/c 2 < m(A+n-) < 2465MeV/c 2 
is satisfied (veto cuts). These two requirements keep 98% 
of the signal, while they remove 85% of B~ -» <_7°(2455)p 
events. The remaining 15% of the background events do 
not peak in the signal AE — tties region. 

The continuum background is reduced by a require- 
ment on the thrust value of the event T < 0.75, where 
we include both charged particles and photons in this 
calculation. The thrust is defined as 



T 



EilPil 



(1) 



where T is the thrust axis defined as the direction which 
maximizes the sum of the longitudinal momenta of the 
particles, and Pi the momentum vector of the i-th par- 
ticle in the CMS. This selection keeps 83% of the signal 
but only 25% of the continuum background, as deter- 
mined from MC simulation and continuum data collected 
40 MeV below the T(4S) energy. 

To further reduce the background from continuum and 
BB events, mainly coming from 77 combinations of low- 
energy, only one B° candidate per event is selected. In 
events with more than one candidate (about 10% of the 
events), first the candidate (s) with the invariant mass 
771(77) closest to the ir° nominal mass are selected. For 
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FIG. 2: Fitted toes distribution without efficiency correction 
(data points); the result of the fit (solid line) and the back- 
ground estimate (dashed line) is shown. 



events with multiple candidates containing the same tt , 
the candidate with the p K~ ir + mass closest to the nom- 
inal A c mass is retained. If there are still multiple B can- 
didates, the candidate with the highest probability of the 
kinematic vertex fit is used. Figure 1 shows a comparison 
between the AE distribution of candidates reconstructed 
in data and in signal MC events, in which signal events 
are obtained by a fit to the tttes distribution in every AE 
bin, as described below. 

The number of reconstructed signal candidates is de- 
termined from a binned x 2 fit to the observed t71es dis- 
tribution shown in Fig. 2. The sum of two Gaussian 
distributions with different means is used to describe 
the signal. The parameters of the two Gaussians are 
fixed to the values obtained from a fit to signal MC 
events. The background is described by the function 
[12] f bg = tt x m ES ^/l-(m ES /m )2x e^M^s/"*) 3 ), 
where mo = 5.289 GeV/c 2 is the kinematic end-point 
value, c a shape parameter left free in the fit, and n is the 
normalization. There are 273 ± 23 signal candidates seen 
in data and the significance of this observation is more 
than lOcr. 

The number of produced signal events used to mea- 
sure the branching fraction is determined by a fit to 
the efficiency-corrected tttes distribution using the same 
parametrization as before. The events are weighted with 
the inverse of the efficiency as a function of the invariant 
mass T7r(/l+7r ). To compute the efficiency the signal MC 
sample is divided in 10 intervals o{m(A+Tr°). For each in- 
terval the tttes distribution is fitted to extract the signal 
MC yield. The efficiency for each interval is computed 
dividing the yield by the number of events generated in 
this interval. The resulting efficiency distribution is then 
fitted by a 4th order polynomial. The averaged signal 
efficiency is 6.0%. 

The weighted data tttes distribution is shown in Fig. 3 
and the fit found 4528 ± 403 signal events (iV S i gna i). The 
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FIG. 3: Efficiency-corrected rriES distribution for B° — > 
A^p-ir (data points). The result of the fit (solid line) and 
the background estimate (dashed line) is shown. 



branching fraction is then calculated as 



750 



B(B 



• 4W; 
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signal 



B(At -> pK-iT 
= (1.94 ±0.17) x 10 



2N 

-4 



B°B° 



(2) 



where the uncertainty is statistical only from the fit, and 
B(A+ -)■ P K-ir+) = (0.050 ± 0.013) [1]. The quantity 
N BOE o = (233.6 ± 2.6) x 10 6 is the number of B°B° pairs 
and B{T(4S) -> B°B°) = 0.5 is assumed. 

To check for peaking background from other B decays 
and random 77 combinations, the analysis is repeated for 
selected samples without mass constraints on the 7r° and 
/1+ mass. The signal yields after subtraction of back- 
ground obtained from the invariant mass distributions of 
the 7T° and A£ are found to be consistent with the default 
analysis. 

The systematic uncertainties are mainly derived from 
studies of data control samples and by comparison of 
data and MC events. The main systematic uncertainty 
arise from differences between data and MC events in the 
AE distribution seen in Fig. 1. The difference between 
the cut efficiency in MC and data, relative to the MC 
one, is used as the systematic uncertainty (4.6%). Other 
systematic uncertainties arise from the veto cuts (3.4%), 
the 7T° reconstruction efficiency (3.0%), the particle iden- 
tification (1.2%), the number of B°B° pairs (1.1%) and 
the reconstruction efficiency of charged tracks (0.9%). To 
determine the uncertainty from the MC model we use to 
generate signal events, these signal events are reweighted 
depending on m(pir°) and a new efficiency function is cal- 
culated. The data toes distribution is then corrected for 
reconstruction efficiencies with this function and fitted 
as before. The difference in the number of signal events 
we use as the systematic uncertainty of the specific MC 
model (2.2%). The systematic uncertainty due to the fit 
is determined by changing the cut-off value of the back- 
ground function by 1 MeV/c 2 (0.50%). The individual 
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FIG. 4: Efficiency corrected distribution of the invariant mass 
m(A^p); points are signal data events; histogram shows signal 
MC events assuming phase space distribution normalized to 
the number of data events. 
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FIG. 5: Distribution of the invariant mass of the A^tt sys- 
tem in the region where the 17^(2455) resonance is expected; 
points are for data with ttibs > 5.272 GeV/c 2 , the curve shows 
the fit. 



contributions to the systematic uncertainty are added in 
quadrature, resulting in the total of 7.1%. 

In Fig. 4, the measured m(A^p) distribution is com- 
pared with a MC simulated one, generated with a phase 
space distribution for the decay to A^pir and normalized 
to the number of data events. To extract the signal distri- 
bution events, the Toes distribution is fitted in every bin 
of m(A^p). There is a clear difference in shape between 
data and simulation, with a clear enhancement at low 
mass, with a significance of 5<r for the first bin, assum- 
ing Gaussian statistics. Such an enhancement is seen in 
many other baryonic B decays and also in baryon produc- 
tion, such as e + e~ — > jAA [13], which proceeds through 
different short-distance processes. 

In Fig. 5, the invariant mass of the A~£ir° combination 



is shown, fitted by a Gaussian function for a possible 
£^~ (2455) signal and by the function n x (m(/l+7r°) — 
[ra(/l+) + to(7t°)]) c to describe the non-resonant frac- 
tion of the signal and background using a likelihood fit. 
The shape parameters for the Gaussian are fixed to the 
parameters obtained from simulated events. The fit re- 
turns N s +, 2455 -< — 3 ± 3 signal events. Therefore, there 
is no evidence for B° — > Z"+(2455)p. The reconstruc- 
tion efficiency for B° ->• £+(2455)p is (1.70 ± 0.05)%. 
Integrating the likelihood function of the fit parameter 



N 



£?(2455) 



> 0, we obtain a Bayesian upper limit at 90% 

confidence level (C.L.) of B(B° -> S+(2455)p)xB(A+ -> 
K-pir+) < 1.5 x 10~ 6 . 

In conclusion, we have observed the decay B° — > 
A+pir and measured the branching fraction as: 

B(B° -> A^pir ) = (1.94±0.17±0.14±0.50)xl0~ 4 , (3) 

where the uncertainties are statistical, systematic, and 
from the A+ branching fraction, A+ — > pK~ir + . The 
ratio of the partial decay width measured here to the 
BABAR measurement of the decay B~ — ► A+pir~ [9] is 

B ^^x^ = 0.61 ±0.09, (4) 



B(B~ -t ATpn- 



1 B° 



where tb- and Tg are the lifetimes of the B mesons. 
This ratio is consistent with the isospin expectation of 
2/3. Given that we don't have evidence for a B° — > 
S^p contribution, wc also compare our B° — > A^pn 
measurement with only the non-resonant contribution to 
the B~ — > A+pir~ decay. We find 



B(B° -> A+pt: ) t b - 

y _Y x -^-=0.80 ±0.11 , (5) 

^ A-c P7T Jnonrcsonant ^B 



B(B 



which is also consistent with the isospin expectation of 
2/3. 

For the resonant subchannel we calculate a 90% upper 
limit of 

B(B° -+ r+(2455)|>) x B(A+ -+ R-pir) < 1.5 x 10" 6 . 

(6) 
The 90% C.L. Bayesian upper limit for the ratio of the 
branching fractions B(B° — > I7+(2455)p) and B~ — > 
r°(2455)p [9] is 



B(B" -j- j;+(2455)p) t_b_ 
B(B- -> E°(2455)p) X ^ 



<0.73 



(7) 



which wc compute by integrating the likelihood profile for 
the ratio of branching fractions over the positive range. 
It is also consistent with the isospin expectation of 2/3. 
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